Astronomy & Astrophysics manuscript no. 
(will be inserted by hand later) 



Line profile variability in the spectra of Oef stars: 
II. HD 1 92281 , HD 1 4442 and HD 1 4434"^ 

M. De Becker and G. Rauw** 

Institut d' Astrophysique et de Geophysique, Universite de Liege, 17, AUee du 6 Aout, B5c, B-4000 Sart Tilman, Belgium 
Received date / Accepted date 

Abstract. We present the very first analysis of the spectroscopic variability of the three rapidly rotating Oef stars HD 192281 
(05(ef)), HD 14442 (05.5ef) and HD 14434 (06.5(ef)). Radial velocities of the He II A 454 1 Une reveal no evidence of binarity 
on time scales of a few days, or from one year to the next, for any of the targets. The He II A 4686 double-peaked emission 
and, to some extent, the H/3 absorption line display significant profile variability in the spectra of all three stars. Data gathered 
during different observing runs spread over six years reveal a rather stable time scale for HD 192281 and HD 14442, whereas 
the variability pattern changes significantly from one year to the other The case of HD 14434 is less clear as no obvious time 
scale emerges from our analysis. In a temptative way to interpret this variability, stellar rotation remains a possible clock for 
HD 192281 and HD 14442. However, currently available models addressing stellar rotation fail to explain some crucial aspects 
of the observed variability behaviour, which appear to be even more complex in the case of HD 14434. 
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1. Introduction 

' Monitoring of the spectroscopic variability has been used ex- 
tensively to study the structure of the stellar winds of hot lumi- 
nous stars (e.g. Fullerton et al. ll9961 Massa et al. ll9951 Kaper 
. et al. ll997l I1999I I. Some objects were found to display cycli- 
■ cal variations in their stellar wind lines on time scales of a 
. few days, possibly related to the rotational period of the un- 
' derlying star. Large scale co-rotating structures in the stellar 
winds due to photospheric features (magnetic fields, non-radial 
pulsations,...) could possibly account for this variability (e.g. 
Kaper et al. l 19991 Rauw et al. l2001> . In order to quantify the 
importance of different processes that could contribute to this 
variability, more stars should be investigated. The objects of 
the scarce category of so-called Oef stars (Conti & Leep |1974l l 
appear a priori as promising targets for such an investigation 
since the double peaked morphology of the He II A 4686 emis- 
sion Une has been attributed to the effect of fast rotation. 
The best studied Oef stars are CPup (04ef) and ACep (06ef). 
The two stars display sttong spectral variability partly at- 
tributed to the effect of rotation (Moffat & Michaud 1981 
Howarth et al.'T995'; Kaper et aI.'T9"99l. In the case of (Pup, 
Eversberg et al. ( 1998 ) further reported stochastic modulations 
of the He II A 4686 line explained by clumps distributed across 
the wind. These authors also confirmed the existence of spec- 
Send offprint requests to: M. De Becker 
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tral modulations on a time scale compatible with the rotation 
period, near the core of the He II A 4686 line. 
Recently, Rauw et al. (2003 Paper I) analysed the com- 
plex behaviour of the He II A 4686 line in the spectrum of 
BD H-60° 2522 (06.5ef). They suggested that rotation was not 
the only clock ruling this variability, and that possibly non- 
radial pulsations were responsible for some short term variabil- 
ity. In this paper, we report the results of an intensive study 
of three other members of the small group of Oef stars, i.e 
HD 192281, HD 14442 and HD 14434. 

HD 192281 (V = 7.6) is an 05(ef) star of the Cyg OB8 
association. On the basis of its large radial velocity quoted in 
the literature, this star was considered as a runaway candidate 
by Gies & Bolton (1986), but the results of their study did not 
confirm this status. The radial velocity of HD 192281 was also 
studied by Barannikov ( 1993) who argued that it was a binary 
with a 5.48 day orbital period, but this result has, so far, never 
been confirmed. 

HD 14442 (V = 9.2) and HD 14434 (V = 8.5) are respectively 
05.5ef and 06.5(ef) stars which were believed to belong 
to the Per OBI association. On the one hand, HD 14442 
has been classified as nitrogen enriched (ON) by Bisiacchi 
et al. ('1982V On the other hand, HD 14434 may have a 
normal nitrogen spectrum (Schild & Berthet 1986). Both stars 
have been reported as blue straggler candidates by Kendall 
et al. dl 9951 1 who considered they could have followed a 
quasi-homogeneous evolutionary track. However, in a subse- 
quent study, Kendall et al. ([1996) found no CNO abundance 
anomalies in the spectrum of HD 14442 and HD 14434. These 
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Table 1. Summary of some properties of HD 192281, 
HD 14442 and HD 14434. References are the following: (1) 
Conti & Leep ( fT974l . (2) Walbo rn ([T972t . (3) Walborn (IT97^ . 
(4) thi s study, (5) Bieging et al. (TWi, (6) Conti & Burnichon 
(QlTSj, (7) Leith erer (1988 1. (8) K endall et al. (1995 ) . (9) Gies 
& Bolt on ( fT986l l. (10) Humphreys (fT978l l. (11) Conti & Ebbets 
(fT977l . 





HD 192281 


HD 14442 


HD 14434 


Sp.T. 


05.5(ef)" 


06ef' 


06.5(ef)'''' 




05Vn((f))^ 


05n(f)p-^ 


05.5Vn((f))p2 




05(ef)'* 


05.5ef 




log(Toff) 


4.67^ 


4.62^ 


4.60'' 


M* (Mo) 


IQ' 


44.548,^ 


44.528." 




93^ 






R(R©) 


15^ 


7.8-16.6**''' 


10.0-12.6*''' 




195 






d (kpc) 


1.78^ 


2.88* 


3.00* 


Association 


CygOB8"' 


PerOBli°'= 


PerOBli°''= 


Vrot sini (kms~^) 


270" 


273" 


400" 



The two values quoted are respectively the spectroscopic and the 
evolutionary masses. 
^ The first value of the radius is calculated for an assumed value 
of the distance modulus, and the second one is derived from the 
spectral type. 

The membership to this association is uncertain (see Kendall et 
al.[T995) 



authors therefore rejected the blue straggler status of these 
stars and questioned their membership in Per OBI, suggesting 
that their distance could be larger than that of this association. 
The three stars have a large projected rotational velocity: 
Conti & Ebbets ( I1977> reported Vrot sini values of 270, 273 
and 400 kms"! for HD 192281, HD 14442 and HD 14434 
respectively. Very few studies have been performed on these 
three stars. In Tabled we summarize the information regarding 
their fundamental parameters that we have gathered from the 
literature. Different spectral types have been assigned to these 
stars and some of them are given in Table[l]for information. 
We shall briefly return to this point in Sect.|3l 

This paper is the second (and last) one of a series devoted to 
the line profile variability of Oef stars. It describes the first line 
profile variability study performed on HD 192281, HD 14442 
and HD 14434. In Sect.|2l, we present the observations and the 
data reduction. Section|3l discusses the mean spectrum and the 
spectral type of our three target stars. Radial velocities and 
equivalent widths are discussed in Sect.|4] Section|5l describes 
the results of our line profile variability analysis, and Sect.|6l 
consists in a discussion of our results and in an attempt of a 
consistent interpretation. The conclusions and prospects for fu- 
ture work are finally given in Sect.E] 

2. Observations and data reduction 

Spectroscopic observations of our three targets were collected 
during seven observing campaigns between September 1998 



and October 2003 at the Observatoire de Haute-Provence 
(OHP). All observations were carried out with the Aurelie 
spectrograph fed by the 1.52 m telescope (Gillet et al. I1994> . 
For the 1998 and 1999 missions, Aurelie was equipped with a 
Thomson TH7832 linear array with a pixel size of 13 /im. From 
2000 on, the detector was replaced by a 2048 x 1024 CCD EEV 
42-20#3, with a pixel size of 13.5 /im squared. 
An overview of the different campaigns is displayed in Table|2| 
Data reduction followed the same approach as described in 
Paper I. Table|2l also specifies for each campaign which star 
was observed, as well as the spectral range covered during the 
run. In the following, we will merge the data collected during 
the two 1999 campaigns (July and August) into one single set 
which will be called Summer 1999, because of the small num- 
ber of spectra taken during these two individual runs. 

3. I\1ean spectrum 

3.1. HD 192281 

The mean blue spectrum of HD 192281 is shown in Fig. [2 The 
data of this spectrum are from the 2001 mission (see Table|2j. 
This part of the spectrum is dominated by absorptions of He I A 
4471, He II AA 4541, 4686 and by the strong H/3 fine. Mg II A 
4481 is clearly present, partially blended with He I A 4471. In 
emission, the most striking feature is due to N III A A 4634-41 
and, to some extent, to C III A 4650 and probably C IV A 4662. 
The two emission lines at AA 4487,4506, identified as S IV lines 
(Werner & Rauch 2001 ), are also weakly present. We empha- 
size the probable presence of N V absorption lines at 4605 and 
4620 A . The presence of these lines in an 05 star is some- 
what surprising. Indeed, these lines are commonly found in the 
spectra of very early (03) stars (Walborn & Fitzpatrick 2000t. 
However, the lines are observed in the spectrum of HD 14947 
(05If) as reported by UnderhiU et al. (1989 1, but not in the 
spectrum of HD15629 (05V) discussed by the same authors. 
Some Si IV emission lines at 4631 A and 4654 A are also 
present. 

The double-peaked structure of the He II A 4686 emission jus- 
tifies the spectral type given by Conti & Leep ( 1974i. Such a 
line shape is supposed to result from a rotating wind (see e.g. 
Petrenz & Puis 1996 1. Using the ratio of the equivalent widths 
(EWs) of the He I A 447 1 and He II 4541 lines (log W = -0.47 
± 0.03) we assign an 05 spectral type, according to the clas- 
sification criterion given by Mathys ( 1988). In the following, 
we will adopt an 05(ef) spectral type for this star, where (ef) 
expresses the prominence of the central absorption component 
in the He II A 4686 fine. 

3.2. HD 14442 

Spectra of HD 14442 have been taken during observing runs 
in July 1999, August 1999, September 2000, September 2001, 
September 2002 and October 2003. The middle spectrum in 
Fig-E shows the same general features as HD 192281 except 
for the N V lines. The central absorption of the He II A 4686 
line appears less pronounced than in the case of HD 19228L 
The strong emission wings observed for HD 14442 confirm the 
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Table 2. Observing campaigns for the study of HD 192281, HD 14442 and HD 14434. The last three columns specify the number 
of spectra obtained for each star during a specific observing run. 



Campaign 


Detector 


Spectral range 


Recipr. disp. 


Res. power 


HD 192281 


HD 14442 


HD 14434 






(A) 


(Amm^^) 










September 1998 


TH7832 


4455 - 4890 


16 


8000 


6 








July 1999 


TH7832 


4100 - 4950 


33 


4000 


7 


5 





August 1999 


TH7832 


4100 - 4950 


33 


4000 


4 


7 





September 2000 


CCD 


4455 - 4905 


16 


8000 





9 





September 2001 


CCD 


4455 - 4905 


16 


8000 


16 


12 





September 2002 


CCD 


4455 - 4905 


16 


8000 


10 


9 


11 


October 2003 


CCD 


4455 - 4905 


16 


8000 


7 


6 


19 




\ I I I 1 \ \ I I 1 I \ [ I I I \ I \ L 

4500 4600 4700 4800 

Wavelength (Angstrom) 



Fig. 1. Mean spectra of our program stars between 4465 A and 4890 A . From top to bottom: HD 1 9228 1 : calculated from spectra 
of the September 2001 mission, between HJD 2452 163.43 and HJD 2452 170.41; HD 14442: calculated from the spectra of 
the September 2002 mission, between HJD 2452518.54 and HJD 2452533.48; HD 14434: calculated from the spectra of the 
October 2003 mission, between HJD 2452916.49 and HJD 2452934.38. 



Oef classification assigned by Conti & Leep (I1974> . 
The equivalent width determined for the He I A 447 1 and He II 
4541 lines allowed us to apply the log(M^') criterion as defined 
by Mathysl 19881. The result (-0.41 ± 0.03) implies an 05.5 
spectral type. In the following of this paper, HD 14442 will be 
considered as an 05.5ef star. 



3.3. HD 14434 

This star was observed during the September 2002 and October 
2003 observing runs. The bottom spectrum of Fig.Q displays 
the same general features as the spectrum of HD 14442, but 
the lines are even further broadened by the fast stellar rota- 
tion. The fact that the N III A A 4634 - 41 emission lines appear 
also significantly broader could imply that these lines are in 
fact formed at the level of the photosphere (rather than inside 



4 



M. De Becker & G. Rauw: Line profile variability in the spectra of Oef stars: part II. 



the stellar wind, see e.g. Mihalas '1973V We further note that 
the fact that N III A A 4634 - 41 lines appear in emission is not 



Radial Velocity (km s ) Hell A 4541 



in agreement with the results of Kendall et al. ( I1996> who re- 
ported the former weakly in absorption, as well as the absence 
of the latter. The double-peaked structure of the He II A 4686 
line dominated by the central absorption, typical of 0(ef), stars 
is obvious. 

From the equivalent width of the He I A 4471 and He II 4541 
lines, we infer log(iy) = -0.19 ± 0.04, corresponding to an 
06.5 type. In the following of this study, we will thus adopt the 
06.5(ef) spectral type for HD 14434, in excellent agreement 



with the spectral type given by Conti & Leep J1974> . 

4. Radial velocities and equivalent widths 

4.1. HD 192281 

As a first step, we have measured the radial velocity (RV) of 
the strongest absorption lines in the blue range. The lines of 
interest here are Hell A 4541 and H/3. Since the Hell A 4541 
line forms probably deep inside the photosphere, it is likely 
to provide the most reliable indicator of the actual stellar RV. 
As shown in Fig.|2] the dispersion in radial velocity reaches a 
maximum of about 10 km s^^, whereas the typical error' in the 
determination of RVs with our gaussian fitting method is about 
8 kms~^. For all our data sets, we find a mean radial velocity 
of —23.1 kms"'^. This leads to a value of about 5.0 kms^'^ for 
the peculiar velocity (after application of the LSR and Galactic 
rotation corrections). This value strongly argues against a run- 
away status for this star Indeed, according to the criteria of 
Stone (1991), the peculiar velocity threshold for a star to be 
a runaway is about 40 kms"'^. It has been suggested by Gies 
& Bolton ( 1986 1 that the runaway status was erroneously at- 
tributed to stars such as HD 192281 because of thekhuge wind 
velocity that yield a large negative RV for those spectral lines 
that are at least partly formed in the wind. The choice of the 
line for RV estimation is thus particularly crucial. 
Gies & Bolton ([1986) measured also the RV of the He II A 4686 
and hydrogen lines for this star. After correction for the veloc- 
ity component due to the solar motion in the LSR frame, they 
found a mean of —26.8 kms~^ for the Hell A 4686 line, and 
a mean of —44.7 kms^^ for hydrogen lines. For comparison, 
we find LSR-corrected RV means of —33.3 and —46.1 kms^^ 
respectively for He II A 4686 and H/3. The reasonable agree- 
ment between our values and those of Gies & Bolton suggests 
that there are no large long-term variations of the RVs. The 
difference of about 18 kms~^ between the mean RV of these 
two lines (about 13 km s~^ in our case) is again explained by a 
difference in the depth of the line formation region within the 
expanding atmosphere. Indeed, the He II A 4686 line is most 
probably formed deeper in the wind, where the temperature is 
higher, and the wind velocity is lower. 

Because of its profile variability, which will be discussed 
hereafter, the He II A 4686 line yields less reliable RVs. Indeed, 
the radial velocities spread over more than 35 km s^^. We did 

' The error on the radial velocity estimate corresponds to tlie 
standard deviation determined for the radial velocity of a Diffuse 
Interstellar Band (DIB) at about 4762 A . 



-35 



-30 
- -35 



2525 2530 
JD - 3450000 



2535 



-40 



Fig. 2. Radial velocities of HD 192281 determined from the 
Hell A 4541 line for the September 2002 observing run. The 
central wavelength was determined by fitting a gaussian to the 
line. The error bar gives an idea of the standard deviation ob- 
tained for the RV data set. 



however not find any systematic trend for this variability in ra- 
dial velocities. 



Barannikov ( I1993> performed a RV variability analysis, the 
only known, for this star. He claimed to find RV variation 
with a semi-amplitude of 16.1 ± 2.4 kms^^ on a 5.48 days 
period, compatible with the existence of a low mass compan- 
ion. We have folded our RV data with Barannikov's 5.48 days 
period but no systematic trend appeared. A Fourier analysis 
was also performed on our RV measurements of the He II A 
4541 line. The resulting periodogram does not show any peak 
with an amplitude larger than 2.4 kms~^. The highest peak 
appears at a frequency of 0.105 d^^, corresponding to a 9.57 
d period. This value is quite close to the period (9.59 d) ob- 



tained by Barannikov (I1993> for his photometric data. In fact, 
Barannikov reported also apparently periodic light variations 
with a peak to peak amplitude of about 0.04 mag. This coinci- 
dence is somewhat puzzling. Although this could be regarded 
as evidence for binarity, we do not consider this to be a plau- 
sible scenario. In fact, the amplitude of the RV variation (^^2.5 
kms^^) at the 9.57 d 'period' is of the order of the uncertainty 
on individual data points. Moreover, it seems extremely dif- 
ficult to account simultaneously for a light curve of amplitude 
^^0.04 mag (peak to peak) and a RV variation inferior to about 5 
kms^^ (peak to peak) within a binary scenario. Assuming that 
the light variations are due to total eclipses of a companion, we 
could be dealing with a secondary of spectral type ~ B3V and 
such an object would produce a RV curve with a peak to peak 
amplitude of > 80 kms~^ (assuming i > 60°). On the other 
hand, to account for the amplitude of the RV curve, we would 
need a companion of mass inferior to about 1 M0 and such a 
star would be far too faint to produce a variation of 0.04 mag 
in the light curve, even if it were still on a pre-main sequence 
track. 

The equivalent width was measured for the He I A 4471, He II 
A 4541 and H/3 lines. The results, as well as the RVs, are col- 
lected in Table|3] Mean values are given for each observing run 
with an estimate of the standard deviation. We see that no sig- 
nificant variability appears from one run to the other Moreover, 
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Table 3. Radial velocities and equivalent widths of the strongest absorption lines in the spectrum of HD 192281. The RV of the 
He I A 4471 line is not quoted because its shape deviates significantly from a Gaussian, leading to a poor estimate of its central 
position. The complex shape of the He II A 4686 line, and its strong variability, prevented us to derive any useful EW for this line. 
The error bars represent the ± 1 cr standard deviation. An estimate of the mean signal-to-noise ratio of an individual spectrum is 
given in the last column. 



Mission 


He I A 4471 




Hen A 4541 


He II A 4686 




H/3 


S/N 




EW 


RV 




EW 


RV 


RV 




EW 






(A) 


(km 


s 


') 


(A) 


(kms-^) 


(km s~ 


') 


(A) 




September 1998 


0.31 ±0.04 


-22.4 


± 


3.0 


0.84 ± 0.04 


-48.9 ± 7.4 








380 


Summer 1999 


0.28 ± 0.01 


-21.5 


± 


2.4 


0.81 ± 0.02 


-43.7 ± 9.2 


-62.8 ± 


3.5 


1.97 ±0.03 


290 


September 2001 


0.31 ±0.01 


-22.9 


± 


2.8 


0.92 ± 0.01 


^4.3 ± 5.0 


-56.5 ± 


2.4 


2.05 ± 0.03 


450 


September 2002 


0.30 ±0.01 


-25.6 


± 


2.9 


0.94 ± 0.01 


^8.6 ± 6.0 


-58.1 ± 


3.9 


2.05 ± 0.02 


440 


October 2003 


0.31 ±0.01 


-13.8 


± 


5.2 


0.88 ± 0.02 


-^1.8 ± 12.8 


-58.1 ± 


3.2 


2.06 ± 0.04 


380 



Table 4. Same as Tabled but for HD 14442. 



Mission 


Hex A 4471 


Hell A 4541 


He II A 4686 






S/N 




EW 


RV 


EW 


RV 


RV 


EW 






(A) 


(kms-i) 


(A) 


(kms-i) 


(kms^^) 


(A) 




Summer 1999 


0.33 ± 0.02 


-30.3 ± 8.2 


0.84 ± 0.02 


-51.5 ± 18.4 


-67.7 ± 3.7 


1.84 ± 0.04 


160 


September 2000 


0.32 ± 0.02 


-32.2 ± 6.7 


0.84 ± 0.02 


-59.7 ± 14.9 


-66.6 ±5.1 


1.91 ±0.03 


340 


September 2001 


0.34 ± O.OI 


-31.8 ± 5.6 


0.84 ± 0.01 


-57.1 ± 23.7 


-64.5 ± 3.8 


1.86 ±0.03 


300 


September 2002 


0.32 ± O.OI 


-34.3 ± 2.9 


0.83 ± 0.02 


-66.5 ± 16.2 


-65.7 ± 4.3 


1.84 ±0.05 


430 


October 2003 


0.32 ± O.OI 


-32.0 ± 10.7 


0.80 ± 0.01 


-57.4 ± 19.6 


-64.5 ± 14.5 


1.84 ±0.05 


330 



Table 5. Same as Table|3] but for HD 14434. 



Mission 


He I A 4471 


Hen A 4541 


He II A 4686 




H/3 


S/N 




EW 


RV 


EW 


RV 


RV 




EW 






(A) 


(kms"^) 


(A) 


(kms~^) 


(km s^ 


') 


(A) 




September 2002 


0.54 ± 0.03 


-15.3 ± 5.7 


0.81 ± 0.03 


^0.2 ± 18.2 


-39.4 ± 


3.5 


2.09 ± 0.03 


530 


October 2003 


0.52 ± 0.02 


-3.9 ± 12.0 


0.84 ± 0.03 


-28.4 ± 13.4 


-35.3 ± 


5.3 


2.10 ±0.05 


320 



no significant variability was found during individual observ- 
ing campaigns. 

4.2. HD 14442 

Gaussians were fitted to the He II A 4541 line to determine the 
radial velocity. The results are illustrated in Fig.lsjwhich shows 
the RV versus heliocentric Julian day for the September 2002 
observing run. We find, for this line, a mean RV of about —32. 1 
kms~^ for the whole data set, which is somewhat different 
from the values reported by other authors such as Conti et al. 



Radial Velocity (km s ') Hell \ 4541 



l fT977t who give -42.1 kms^^ or Underhill & Gilroy ( 1990) 
who published a value of —40.6 kms^^, for averages of sev- 
eral line measurements. As for the differences in RVs found for 
different lines, the discrepancy between our results and those 
of the authors cited hereabove can probably be explained by 
the velocity gradient in the stellar wind. The maximum differ- 
ence in RV over the five data sets is only about 30 kms^^. A 
Fourier analysis was performed on these RV data but no peak 
with an amplitude larger than 4.4 km s^^ emerged from the pe- 
riodogram. 

RVs were also determined for the He II A 4686 and H/3 lines. 




2525 2530 
JD - 2450000 



Fig. 3. Same as Fig.|2l but for HD 14442. 



For the first one, RV values vary around a mean of —57.0 
kms^^, whilst for H/3, this mean is about —66.3 kms^^. As 
for HD 192281, the mean RV is more negative in the case of 
the hydrogen line. 

The equivalent width was determined for the same lines as 
for HD 192281. Once again, no significant variability is found 
during the runs, or from one observing run to the other. 
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Radial Velocily (kii 



) Hell A 4541 



2525 2530 
JD - 3450000 



Fig. 4. Same as Fig.|2l but for HD 14434. 
4.3. HD 14434 

Radial velocities measured on the same lines as for HD 192281 
and HD 14442 are given in Table|5] We emphasize that these 
values are subject to larger uncertainties due the difficulty to fit 
gaussians to the very broad absorption of this star: the actual 
line profiles are not gaussian and undergo variations which are 
discussed in more details in Sect. l5.3l The rather large values 
of the standard deviations given in Tablels] reflect these diffi- 
culties. The largest error bars are obtained for the He II A 4686 
line, and are undoubtedly due to the strong variability of this 
line. 

Fig.EJshows the radial velocities measured on the He II A 4541 
line during the September 2002 observing run as a function of 
time. The dispersion is rather large, but no systematic trend can 
be detected. Mean radial velocities for the September 2002 and 
October 2003 campaigns overlap within their large error bars. 
The radial velocity reported by Conti et al. (1977 ), i.e. -11.7 
kms^^, is compatible with the value obtained for the He II A 
4541 line. No indication of binarity on time scales of several 
days, or from one year to the next, has thus been revealed by 
our time series. 

Equivalent width are also quoted in Table|5] These EWs are 
stable on the time scales investigated in our study. 



5. Line profile variability 

The first step of our line variability analysis (see also Rauw et 
al. l2003t is the identification of the lines which vary signifi- 
cantly with time. For this purpose, we used the time variance 



spectrum (TVS) as developed by Fullerton et al. (I1996> . We 
evaluated the signal-to-noise ratio of individual spectra over a 
line free region between 4560 A and 4590 A . This estimate is 
necessary for the determination of the 99 % significance level 
of the variability found with the TVS. 

The procedure was applied to each dataset. The lines that 
show significant (at the 99 % confidence level) variability 
are essentially He II A 4686 and H/3. The detailed line profile 
variability of these two lines is described hereafter. 

In order to quantify the variability of the various spec- 
tral lines, we performed a 2D Fourier analysis using the tech- 



nique described by Rauw et al. (12001 > : at each wavelength step, 
the detection of candidate frequencies is performed by an un- 



equally spaced data Fourier analysis as described by Heck et 
al. (1985 1. The main peaks (ly) and their aliases (essentially 
l—i/) are identified and are considered for further analysis. As a 
next step we 'prewhite' the data for the variations at frequency 
I' following a method already described by Rauw et al. (2001) 
and used in Paper I. The frequencies which give the best results, 
i.e. a periodogram devoid of significant peaks after prewhiten- 
ing, are considered as good candidates. Ideally, a physically 
significant frequency should meet the following criteria: 

1 . A stable clock ruling the variability should be found in the 
data of a given line from more than one observing cam- 
paign. 

2. The frequency should not be found in the periodogram of 
a constant line (typically an interstellar feature), otherwise 
such a frequency probably reflects an artificial variability 
affecting the entire spectrum, and should not be associated 
with a physical process. 

3. The frequency should appear, within the error bars, in the 
periodograms of more than one variable line if more than 
one line varies significantly with time. 

These criteria concern in fact only the most ideal cases and 
have to be relaxed in more complex situations (like the case of 
BD H-60° 2522, Paper I). Their strict application would proba- 
bly lead to the rejection of potentially meaningful frequencies. 
The last criterion, in particular, is certainly not mandatory be- 
cause we cannot state that a given variability mechanism will 
affect all lines in a same manner 

5.1. HD 192281 

An overview of the observing runs is given in Table|6l 
Generally, one spectrum was obtained during each night the 
star was observed, except for the 2001 run where at least 2 
spectra were taken each night with a mean time separation of 
about 0.025 d. 

He II A 4686. As shown in Table|6l five data sets were avail- 
able for the study of this line. TVS and Fourier analyses were 
performed on the five time series, as well as on a combination 
of all our data sets. For a better visualization of the behaviour 
of this line, Fig.|5] shows the mean spectrum and the TVS, as 
well as the 2-D power spectrum as a grey scale representation 
for the September 2001 observing run, and in the case of the 
combined data set. The periodogram of the combined data set 
displays some kind of fine structure due to the strong irregular- 
ity of the sampling. The TVS calculated for all our individual 
data sets reveals that the level of the variability changes from 
one observing run to the other The variability affects essen- 
tially two separate parts of the line, respectively the violet and 
the red wings. 

The most striking and recurrent feature in these peri- 
odograms is the presence of aliases at frequencies of about 0.35 
and 0.65 d^^, corresponding respectively to periods of about 
2.86 and 1.54 d. Figure|6l shows the mean periodogram ob- 
tained for the 2001 run between 0.0 and 5.0 d~^, along with 
the residual periodogram obtained after prewhitening for a fre- 
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Table 6. Observing runs used for the line profile variability 
study of HD 192281. The first column gives the name of the 
campaign as used in the text, the second column yields the 
amount of spectra taken during each run, the third specifies 
which line was studied with the corresponding data set, the 
fourth one gives the time elapsed between the first and the last 
spectrum of the run, and finally the last column provides the 
natural width of a peak of the power spectrum taken as 1/AT. 



Mission 


Number of 


Line 


AT 






spectra 




(d) 




Sept. 1998 


6 


He II A 4686 


7.977 


0.125 


Sum. 1999 


11 


He II A 4686 


33.140 


0.030 






H/3 






Sept. 2001 


16 


He II A 4686 


6.979 


0.143 






H/3 






Sept. 2002 


10 


He II A 4686 


14.984 


0.067 






H/3 






Oct. 2003 


7 


He II A 4686 


15.910 


0.063 



H/3 



See additional GIF figures 



Fig. 5. Upper panels: mean spectrum and square root of the 
time variance spectrum of the He II A 4686 line of HD 192281 
between 4675 A and 4695 A respectively for the September 
2001 observing run and for the combined data set. The 99 % 
significance level of the TVS^^^ is indicated by the horizontal 
dotted line. Velocities are expressed in kms~^. Lower panels: 
the corresponding two-dimensional power spectra. The dark- 
est regions stand for the highest peak intensities of the peri- 
odogram. Frequencies and wavelengths are expressed respec- 
tively in d^^ and in A (see figure 0548f5.gif). 



quency of 0.64 d^^ between 4687.5 and 4688.5 A . The spectral 
window is displayed in the bottom panel. We emphasize that 
the time analysis performed on such a narrow wavelength band 
does not provide any result representative of the behaviour of 




1 2 3 4 5 



Fig. 6. Upper panel: mean power spectrum obtained for the 
He II A 4686 line of HD 192281 between 4687.5 A and 4688.5 
A for the September 2001 observing run, up to 5.0 d~^. Middle 
panel: residual power spectrum obtained after prewhitening 
with a frequency of 0.64 d^^. Bottom panel: power spectral 
window showing secondary peaks which explain the asymmet- 
rical shape of the peaks of the power spectrum. 



the whole line, but it allows a clear identification of the time 
scale of the variability observed in the part of the line where the 
TVS reaches its highest level. Prewhitening shows that most of 
the power of the periodogram can be accounted for by aliasing 
of this 0.64 d^^ frequency (or its 1 — alias). 

The periodogram of a Diffuse Interstellar Band (DIB) at 
4762 A has been calculated for each data set, and only the 
September 1998 campaign gives peaks situated near these can- 
didate frequencies, but with much lower intensities. 
The variability pattern of the He II A 4686 line changes with 
epoch. Indeed, the strongest variability and hence the v\ and 
\~v\ peaks appear essentially in the red and the violet (with 
lower intensity) wings in the September 1998 and September 
2001 periodograms, but appear only in the violet wing in 
Summer 1999. The width of the wavelength range displaying a 
significant variability with frequency v\ or X—vi changes also 
from one epoch to the other. A summary of all the frequencies 
is given in Table|9](Sect.|6}. 

H/3. The same approach as for the Hell A 4686 line was 
adopted for H/3. Our data allow us to perform this line pro- 
file variability analysis only for the 1999, 2001, 2002 and 2003 
campaigns. In this case, the TVS and the 2-D power spectra re- 
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Table 7. Same as Table|6j but for the observing runs used for 
the Hne profile variability study of HD 14442. 



Mission 


Number of 


Line 


AT 


Afnat 




spectra 




(d) 




Sum. 1999 


11 


He II A 4686 


34.039 


0.029 






H/3 






Sept. 2000 


9 


He II A 4686 


10.967 


0.091 






H/3 






Sept. 2001 


12 


He II A 4686 


7.030 


0.142 






H/3 






Sept. 2002 


9 


He II A 4686 


14.942 


0.067 






H/3 






Oct. 2003 


6 


He II A 4686 


17.957 


0.056 



H/3 



veal that the variability is concentrated in the core of the line. 
The Summer 1999 data set analysis leads once again to fre- 
quencies of about 0.33 or 0.67 d^^. This observation raises thus 
the degree of confidence concerning these frequencies akeady 
noted for the He II A 4686 line. However, the situation is less 
clear for the September 2001 periodogram. Indeed, the peak 
frequencies found in this case are near 0.14 and 0.46 d~^. For 
the September 2002 campaign, the variability appears shifted 
to the red wing of the line, and we are not able to definitely 
distinguish the peaks from those which appear in the DIB peri- 
odogram. 

He I A 471 3. Even if no variability is seen for this line until 
2001, the September 2002 observing campaign reveals a sig- 
nificant TVS signal in the red wing. The existence of this vari- 
ability turns out to be robust in the sense that it remains if we 
omit individual spectra one at the time and repeat the analysis 
on the remaining data set. Although the periodogram is dom- 
inated by a peak at ~0.3 d~^, prewhitening at this frequency 
yields no satisfactory result. 

5.2. HD 14442 

Table0 gives some information about the observing runs. 
During the September 2001 mission, two spectra were taken 
each night with a mean separation between the two spectra of 
about 0.025 d. 

He II A 4686. In Fig.|8] we show the results obtained for the 
September 2001 campaign, as well as for the combined data 
set. The analysis of the individual data sets reveals that the 
TVS obtained for the September 2002 campaign shows a more 
complex structure, with an additional varying component in the 
blue emission peak. To check the robustness of this feature, we 
have removed each individual spectrum one at the time from 
our September 2002 time series and performed the variability 
analysis on the remaining data. This additional component in 
the TVS appears to be stable. The situation is even more com- 
plex in the case of the October 2003 data set, where we can 
distinguish four distinct varying parts in the line, but it could be 




1 2 3 4 5 

u id-') 



Fig. 7. Upper panel: mean power spectrum obtained for the 
He II A 4686 line of HD 14442 between 4682 A and 4684 A 
for the September 2001 observing run, up to 5.0 d~^. Middle 
panel: residual power spectrum obtained after prewhitening 
with a frequency of 0.60 d^^. Bottom panel: power spectral 
window showing secondary peaks which explain the asymetri- 
cal shape of the peaks of the power spectrum. 



due to a rather poor sampling. Due to the poor weather condi- 
tions we obtained indeed only 6 spectra during that campaign. 
All our data sets (except perhaps in Summer 1999 where the 
variability is much less significant) yield a relatively stable fre- 
quency between 0.35 and 0.40 d^^ (or 0.65 and 0.60 d"^ for 
their respective aliases) near the line core. These frequencies 
give very good prewhitening results and are not the same as 
those of the main peaks appearing for the DIB at 4762 A . 

Fig0 shows the mean power spectrum obtained for the 
September 2001 mission between 4682 A and 4684 A. The 
middle panel of this figure illustrates the residual periodogram 
obtained after prewhitening with a frequency of 0.60 d~^. 
It appears that the aliasing of this frequency is sufficient to 
account for all the peaks of the power spectrum. A similar 
result is achieved if the prewhitening is performed with a 0.40 
d^^ frequency. 

Other peaks appear also in our periodograms. For example, 
a frequency of 0. 1 1 d^^ (and its aliases) is clearly identified in 
the red wing for the Summer 1999 campaign. Table|9]in Sect.|6l 
gives a more complete overview of the frequencies emerging 
from the analysis of all our time series. 
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See additional GIF figures 



Fig. 8. Same as Figure|5] but for HD 14442 (see figure 
0548f8.gif). 



See additional GIF figures 



Fig. 9. Same as Figure Elbut for the H/3 line between 4845 A 
and 4875 A (see figure 0548f9.gif). 



H/3. No significant variability was found for this line in the 
Summer 1999 data. Remember that this data set has the poorer 
signal-to-noise ratio, and also the poorer sampling due to the 
spread of our observation over 34 days. However, significant 
variations were found for the other observing runs, especially 
for the September 2002 and the October 2003 campaigns where 
it reaches its highest level. The two-dimensional periodograms 
in Fig.|9] show variability patterns which are very similar to 
those obtained for the He II A 4686 line. The simplest case ap- 
pears for the 2001 data set, where the main feature is com- 



pletely reproduced by a 0.39 d^^ (or 0.60 d^^) frequency. 
This is very interesting if we remember that it corresponds 
quite well to the frequencies obtained for the Hell A 4686 
line. The situation is more complex in the case of the 2000 
and 2002 campaigns. After analysis, it appears that the peri- 
odograms are partly explained by some complex long term be- 
haviour probably due to the sampling or normalization irreg- 
ularities. Nevertheless, frequencies of about 0.35 d^^ (or 0.65 
d^^) obviously contribute to the power spectra obtained in both 
cases. The October 2003 periodogram is mainly dominated by 
peaks at about 0.41 d^^ (0.59 d^^). 

He II A 4541 . In the 2002 campaign, the good quality of the 
data allowed us to find a significant (albeit weak) line profile 
variability in the case of the He II A 4541 line, especially near 
the core of the line with some weak contributions in the red 
wing. The analysis of the periodogram revealed that frequen- 
cies between 0.36 and 0.4 d^^ (or their l—i/ aliases) contribute 
significantly to the variability of this line. This range of fre- 
quencies matches very well the frequencies obtained for the 
He II A 4686 and the H/3 lines discussed hereabove. This result 
reinforces our trust in the physical origin of the frequencies 
found for the variations in the spectrum of HD 14442. 

5.3. HD 14434 

HD 14434 was observed during only two observing runs (see 
Table|8}. In order to search for rapid variations, the star was 
intensively observed during the night from October 16 to 17, 
2003: ten spectra spanning 6 hours were obtained. The TVS for 
this specific night computed over the whole spectrum did not 
reveal any significant variability. Consequently, we conclude 
that HD 14434 does not display a variable behaviour on short 
time scales (a few hours) in the spectral range we investigated. 
However, over time scales of several days, a significant vari- 
ability is detected essentially for the He II A 4686 line in both 
data sets, and for H/3 in October 2003. The behaviour of these 
two lines was investigated and the results are described here- 
after The high quality of our September 2002 data also allowed 
us to detect some variability for the He I A 4471 line, but with 
a level which is too low for a Fourier analysis. 

Hell A 4686. In Fig.[lO| the TVS shows that a significant 
variability is detected for the He II A 4686 line in both data 
sets. The variability level seems higher in October 2003 than 
in September 2002. The variations are not restricted to a given 
part of the line, but extend over a rather broad domain. Lower 
panels of the same figure display the two-dimensional power 
spectra obtained through Fourier analysis. The behaviour of the 
line is rather complex, and the periodograms reveal that dif- 
ferent parts of the line apparently display different variability 
trends. A more detailed analysis, i.e. by focusing on individual 
parts of the line, yields different time scales. 
The September 2002 power spectrum reveals that the dominant 
frequency is about 0.91-0.92 d^^. This value corresponds to the 
darkest spots of the lower left panel of Fig.fTO] We note that the 
validity of the corresponding periods (close to one day) can not 
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Table 8. Observing runs used for the line profile variability 
study of HD 14434. 



Mission 


Number of 


Line 


AT 






spectra 




(d) 




Sept. 2002 


11 


He II A 4686 


14.936 


0.067 






H/3 






Oct. 2003 


19 


He II A 4686 


17.891 


0.056 



H/3 



be totally warranted by our sampling. Other peaks contribute 
also to the power but the prewhitening failed to reveal unam- 
bigously which ones were to be selected to identify the related 
variability time scales. In October 2003, the situation is also 
rather complex though the prewhitening indicates that frequen- 
cies at about 0.35-0.4 d^^ dominate the power spectrum. We 
note also that the same analysis performed on the DIBs did not 
reveal significant variability due to normalization eiTors that 
could disturb our variability analysis. 

H/3. This line displays also a variable behaviour, but only at a 
rather low level. The detection of this variability is most signif- 
icant for the October 2003 data set, even though data of better 
quality were obtained during the September 2002 campaign. 
In October 2003, the variable part of the line could be divided 
into a red and a blue side. The red part seems to be ruled by 
a combination of frequencies of about 0.14 and 0.35 d~^ (or 
their 1 — aliases), or by a frequency of about 0.44 d~^. In the 
blue part, the possible frequencies are close to 0.3 or 0.35 d^^ 
(or their 1 — v aliases). Even though the frequencies of about 
0.35 d^^ appearing for the He II A 4686 and Hf3 lines could be 
related to a common physical process, we did not obtain suf- 
ficiently consistent results to give a detailed description of the 
variable behaviour of these lines. 



6. Discussion 

At this stage, it is interesting to recall the different time scales 
which emerged from our line profile analyses for the three stars. 
Tab le|9] provides the frequencies found (as well as the corre- 
sponding periods) for the different lines, and for each data set. 
Because of the ambiguity between aliases (t/ and l—i^), they 
are both given in most cases, with the one of lowest frequency 
in brackets. 

It appears from Table|9]that some variability time scales 
are found in different data sets of a same star. This is obvious 
for HD 192281 and HD 14442. This recurrence in some 
frequencies suggests that a rather stable process is at work in 
the variability we observe. However, the case of HD 14434 
is much less conclusive. Indeed, our data sets did not reveal 
any obvious time scale which could be responsible for a large 
part of the variability detected for the He II A 4686 line, and to 
some extent for H/3. This points to a more complex situation 
than for HD 192281 and HD 14442. Nevertheless, probably 
the most important result arising from our line profile analysis 
of HD 14434 is the absence of variations on time scales of 



See additional GIF figures 



Fig. 10. Upper panels: mean spectrum and square root of the 
time variance spectrum of the He II A 4686 line of HD 14434 
between 4675 A and 4695 A respectively for the September 
2002 and the October 2003 observing runs. The 99 % signif- 
icance level of the TVS^^^ is indicated by the horizontal dot- 
ted line. The differences in significance levels are due to dif- 
ferences in signal-to-noise ratios (see Table|5}. Velocities are 
expressed in kms^^. Lower panels: the corresponding two- 
dimensional power spectrum. The darkest regions stand for the 
highest peak intensities of the periodogram. Frequencies and 
wavelengths are expressed respectively in d^^ and in A (see 
figure 0548fl0.gif). 

several hours. This is especially important in the framework 
of scenarios where non-radial pulsations may play a role (see 
Sect. 16. 2> . The non-detection of short time scale variations 
suggests that pulsations do not have a major impact in this 
Oef star, contrary to what was put forward by the study of 
BDh-60°2522 (Paper I), (Pup (Reid & Howarth [T996> and 
ACep(de Jong et al. 1999 1. 

In addition, there are several points that have to be consid- 
ered in a model for the variability discussed in this study. First, 
the stars are rapid rotators and display a double peaked He II A 
4686 emission line. Second, the 'periodic' variability observed 
during a specific observing run affects only part of the line pro- 
file. Third, the variability pattern changes from one epoch to 
the other and more specifically, the wavelength range display- 
ing the periodic behaviour can shift from one wing of the line 
to the other Finally, the variability level changes as well from 
one oberving run to the next. 

6. 1. Rotational modulation 

Since Oef stars are fast rotators, a straightforward question is 
whether a stable variability time scale as observed in the case 
of HD 192281 and HD 14442 could be linked to stellar rotation. 
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Table 9. Summary of the frequencies emerging from our line profile analyses of HD 192281, HD 14442 and HD 14434. The 
frequencies are given with their l—v aliases in brackets. When two frequencies are quoted the aliases are not specified. When 
more than one value is provided for a line of an observing run, it indicates the results obtained for different wavelength regions 
of the same line, from the blue to the red. ND means that 'no data' were available for that specific line and NV stands for 'no 
significant variability'. 







HD 192281 


HD 14442 


HD 14434 


Observing run 


Line 


V 


P 


V 


P 


V 


P 








(d) 


(d-^) 


(d) 




(d) 


September 1998 


He II A 4686 


63 fO 37^ 


1 59 O 


ND 


ND 


ND 


ND 


Summer 1999 


He II A 4686 






6(S (Ci 'KVt 




ND 


ND 












1.69 (2.44) 








H/3 






NV 


NV 


ND 


ND 


September 2000 


He II A 4686 


ND 


ND 


6? (0 38'! 


1 61 (1 63'! 


ND 


ND 










69 ro 3 n 


1 4s n T\\ 














0.76 (0.24) 


1.32 (4.17) 








H/3 


ND 


ND 


0.65 (0.35) 


1.54 (2.86) 


ND 


ND 


September 2001 


He II A 4686 




1 '^i n 86"! 




1 64 (1 sn~i 


ND 


ND 






64. (Vi ^6"! 


1 "Sfi n 78"! 

1 . JU 1 












H/3 


1 4 * 46 


7 14 & 2 17 


61 rO 39') 


1 64 C2 S«\ 


ND 


ND 


September 2002 


Hen A 4541 


NV 


NV 


0.75 (0.26) 


1.33 (3.85) 


NV 


NV 










0.63 (0.37) 


1.59 (2.70) 








He II A 4686 


0.64 (0.36) 


1.57 (2.78) 


0.63 (0.37) 


1.59 (2.70) 


0.91 


1.10 






0.65 (0.36) 


1.54 (2.78) 


0.66 (0.34) 


1.52 (2.94) 


0.92 


1.09 






0.67 (0.33) 


1.49 (3.03) 


0.77 (0.23) 


1.30 (4.35) 


0.70 & 0.91 


1.43 & 1.10 










0.61 (0.39) 


1.64(2.56) 








H/3 


NV 


NV 


0.65 (0.35) 


1.54 (2.86) 


NV 


NV 










0.77 (0.23) 


1.30 (4.35) 














0.64 (0.36) 


1.56 (2.78) 






October 2003 


He II A 4686 


0.67 (0.33) 


1.49 (3.03) 


0.61 (0.39) 


1.64(2.56) 


0.15 & 0.60 


6.67 & 1.67 






0.59 (041) 


1.69 (2.44) 






0.65 & 0.82 


1.54 & 1.22 














0.66 & 0.83 


1.52 & 1.20 




H/3 


0.13 & 0.26 


7.69 & 3.85 


0.57 (0.41) 


1.75 (2.43) 


0.44 


2.27 



0.14 & 0.74 7.14 & 1.35 
0.34 & 0.74 2.94 & 1.35 



To address this issue, we must first derive some constraints on 
the rotational period of these stars. An upper limit on the rota- 
tional period can be obtained from the observed projected rota- 
tional velocity. 



p _ 27ri?eq 

ma.x — 



Vrot sini 



(1) 



where i?cq is the stellar radius at the equator (see TableO. 
A lower limit can be evaluated as the inverse of the critical 
rotational frequency, beyond which the centrifugal acceleration 
overwhelms the effective gravity acceleration. This condition 
yields the relation 



p . 

mm 



27ri?: 



3/2 



(GAf, (1-r)) 



1/2 



(2) 



where G is the gravitational constant, Af, the mass of the star, 
and r the ratio of the luminosity of the star to the Eddington 



relation given by Lamers & Leitherer ( 11993 



luminosity. This parameter was evaluated with the following 



(3) 



r = 7.6610"^CTo 



Mq 



where all parameters have their usual meaning, with the elec- 
tron scattering coefficient de taking the value 0.34 cm^ g~^. 



6.1.1. HD 192281 and HD 14442 

Let us first consider the case of HD 192281. We adopt 
K-ot sin i =21 Q km s^^ as obtained by Conti & Ebbets J1977> . 
With the values taken from Leitherer ( 1988 1 for the radius (15 
Rq), the luminosity (lO^-^'* L©) and the stellar mass (70 M©), 
we find that the rotation period should be in the range from 
0.98 to 2.81 d. Those boundaries correspond respectively to 
frequencies of 1.02 and 0.36 d~^. 
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We see in Table|9] that most of the periods range between 
1.40 and 1.59 d (with ahases between 2.70 and 3.44 d). It is 
noticeable that most of the peaks of our periodograms lie so 
close to each other, and that the results are similar for several 
epochs. From the limits discussed hereabove, we find that most 
probably only the group of periods near 1.5 d could possibly 
be related to the rotational period. 

The case of HD 14442 is even more interesting. By the use 
of equations Q to (|3jl, we calculated the constraints on the 
rotational period. We adopted a luminosity of log(L/LQ) = 5.62 
intermediate between the two values given by Kendall et al. 
J1995> . For the radius, we chose the value 12.4 Rq corre- 
sponding to the luminosity given hereabove, and the Toff from 
Tabled Kendall et al. ( 1995J estimated the mass of HD 14442 
by combining the stellar radius with atmospheric parameters, 
and by plotting its bolometric magnitude versus Tcff on a 
HR-diagram and interpolating with theoretical evolutionary 
tracks. They found two values of respectively about 44 Mq 
and 54 M0. We use an intermediate value of 49 M0. We find 
limits on the rotational period of Pmin = 0.82 d and Pmax = 
2.29 d. These periods correspond to frequencies of t^max — 
1.22 d-i and z^,ni„ = 0.44 d"!. 

As shown in Tabled the results are quite consistent from one 
line to the other (He II A 4686 and H/3), and from one epoch 
to the other. The first group of periods lies between 1 .52 and 
1.69 d, and the second one between 2.44 and 3.03 d. Since the 
rotational period of HD 14442 should range between 0.82 and 
2.29 d, only the first group of periods could possibly be related 
to a rotational modulation. 
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Fig. 11. Phase coverage of our data sets considering the pro- 
posed periods of 1.5 and 1.6 d respectively for the variability 
of HD 192281 (upper panel) and HD 14442 (lower panel). The 
zero time for phase calculation is arbitrarily chosen at HJD 
2450000 for both stars. Each horizontal series of points cor- 
responds to an observing run whose date is indicated on the 
vertical axis labelled in years. 



riods becomes satisfactory when all our data sets are combined. 



Therefore, assuming that the line profile variability dis- 
cussed here is indeed related to rotation, we tentatively infer 
rotational periods of about 1.5 and 1.6 d for HD 192281 and 
HD 14442 respectively. The uncertainties on these periods 
would be of the order of 0.1 d. From the projected rotational 
velocities and the stellar radii (see Tabled and Sect. 3), we 
can infer inclination angles of 30—35° and 41—47° for the 
rotational axis of HD 192281 and HD 14442 respectively. 
An important issue to address is whether our data sets provide 
a sufficient sampling of the periods that apparently rule the 
line variability in these two stars. Fig^J shows the phase 
coverage of all the data sets used for this study. These plots 
were obtained assuming periods of 1.5 and 1.6 d respectively 
for HD 192281 and HD 14442. In the case of HD 192281, we 
see that most points cluster around three phases. This is due 
to the fact that a period of 1 .5 d cannot be fully sampled from 
a single observatory on the ground, whatever the number of 
spectra obtained. However, as stated in the previous paragraph, 
the confidence interval for this period ranges from 1 .4 to 1 .6 
d. Except around the central value of this domain, the phase 
coverage appears to be satisfactory. The case of HD 14442 
is less problematic. With a period of about 1.6 d, we are 
not confronted to the same difficulty as with HD 192281. 
For each mission, the phase coverage is satisfactory, except 
for values close to the lower boundary of the confidence 
interval. For both stars, the phase coverage for the selected pe- 



Fast rotation had been suggested to produce so-called wind 
compressed zones (Bjorkman & Cassinelli 1993 1 near the stel- 
lar equator. However, more recent models suggest that non ra- 
dial radiative forces prevent the formation of such an equato- 
rial compressed wind (Owocki at al.'1996^ and tend rather to 
focus the wind towards the poles. The non-spherical shape of a 
rapidly rotating star will also lead to enhanced mass loss at the 
pole through gravity darkening (Maeder 19991. Still, an equato- 
rially confined wind could be produced by the effect of a stellar 
magnetic field. 

The rotational broadening of absorption lines in the spectra of 
O stars renders a direct measurement of their magnetic field ex- 
tremely difficult. Nevertheless, dipolar magnetic fields of 360 
G and 1 100 G have been measured for /3 Cep and 9^ On C re- 
spectively using spectropolarimetric techniques (Donati et al. 
^001 2002 1. These magnetic field strengths are sufficient to 
confine the stellar winds into the region of the magnetic equa- 
tor Theoretical models for such confined corotating structures 
have been developed by Babel & Montmerle ( I1997> and very 
recently by ud-Doula (20021. The latter author considered the 
case of aligned magnetic and rotational axes. He found that 
a corotating wind structure forms if the rotational velocity 
reaches a significant fraction of the critical velocity. According 
to ud-Doula (2002 1, the corotating structure in the wind con- 
sists in a density enhancement of the plasma with a lower out- 
flow velocity. This density enhancement tends to increase as the 
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rotational velocity approaches the critical value, and increases 
also as the magnetic field strength increases. In our case, the ro- 
tational velocities reach several tenths of the critical velocities 
and a magnetic field of a few hundred Gauss could be suffi- 
cient to trigger a corotating wind structure. Such a corotating 
structure might explain the double-peaked morphology of the 
He II A 4686 emission in Oef stars. Moreover, if the magnetic 
field is tilted with respect to the rotational axis, we expect the 
corotating confined wind to produce a rotational modulation 
of the spectral lines formed in the wind. Such an 'oblique mag- 
netic rotator' model was put forward e.g. by Moffat & Michaud 
( fT98Tt . Stahl et al. ( IT996t and Rauw et al. (2001) to explain the 
line profile variations of Pup, d^Oii C and HD 192639 respec- 
tively. 

The inclination of an oblique corotating wind structure with 
respect to our line of sight is expected to vary with a period 
equal to the rotational period. This general scenario could lead 
to a large number of variability patterns depending on the in- 
clination angle of the rotation axis on the one hand, and the 
tilt angle between the magnetic and the rotational axes on the 
other hand. The larger the tilt angle, the larger the contrast be- 
tween the observed configuration of the confined corotating 
wind structure at different phases. The radial extent of the coro- 
tating wind would depend on the wind velocity and the mag- 
netic field strength. Therefore, lines formed over a large radius 
would be less affected by the rotational modulation than those 
formed in the innermost part of the wind where the effect would 
be most pronounced. 

The oblique magnetic rotator model might account for the 
recurrence of a stable clock in the line profile variability. 
However, as pointed out above, the variability with this time 
scale only occurs over a limited range of wavelength and shifts 
from the red wing to the blue wing (and vice versa) from one 
year to the next. If our data provide a good sampling of the 
phenomenon, this fact cannot easily be explained by the simple 
tilted magnetic rotator model described hereabove. 

6.1.2. HD 14434 

In the case of HD 14434, stellar rotation does not appear a pri- 
ori as a good candidate to interpret the line profile variations 
that we observed. Our data did not allow us to clearly identify 
frequencies responsible for the variable behaviour investigated 
in Sect. l5.3l However, we should notice that the search for a 
reccurent variabihty time scale requires more observations of 
this star. 

6.2. Other scenarios 

Beside the rotational modulation discussed above, other 
factors are likely to influence the variability of the objects 
studied in this paper Indeed, although the main time scales 
of the variability of HD 192281 and HD 14442 are broadly 
compatible with the rotation period, some aspects do not fit 
this interpretation. This is even more true for the cases of 
HD 14434 and BD H-60° 2522 (Paper I). 



First, in the framework of the simulations of ud-Doula 
( 2002J . it is shown that for large wind confinement parameters, 
i.e. large magnetic field strengths, the density structure near 
the surface becomes quite complex and undergoes a compli- 
cated variability behaviour The occurence of this phenomenon 
at distances from the star corresponding to the formation region 
of the He II line could have a strong effect on its shape. 
Second, we know that the best studied Oef stars (CPup and 
ACep) undergo non-radial pulsations (see Reid & Howarth 
1996 de Jong et al. l999'). Cranmer (1996) considered the 
propagation of non-radial pulsations (NRPs) into accelerated 
stellar winds, and showed how this could lead to large scale 
structures likely able to generate significant line profile modu- 
lations. The study of NRPs in rotating early-type stars is also 
considered by Townsend ( I1997> . According to the latter author, 
rotation compresses pulsational activity towards the equator. 
NRPs and probably multimodal pulsations are believed to trig- 
ger mass loss enhancements (or depletions) possibly responsi- 
ble of the apparition of structures in the stellar wind. Moreover, 
the combined effects of rotation and strong magnetic fields on 
NRPs is not well understood and can lead to complex cases 
like the so called 'off-axis modes' encountered for Ap stars 
for which the axis of pulsation does not coincide with the ro- 
tation axis (Townsend 1 1 997> . These considerations could pro- 
vide the ingredients for an alternative scenario. The interplay 
between NRPs, rotation and the stellar wind were proposed to 
influence the variability behaviour of BD H-60° 2522 (Paper I), 
because of the lack of stability in the variability pattern in that 
case. However, at least in the case of HD 14434, a series of 
spectra obtained with about 30 minutes separation over more 
than 6 hours in October 2003 did not reveal significant short- 
term variability for any of the spectral lines in our wavelength 
domain. This result argues against a major role of pulsations 
at least for HD 14434, but high signal-to-noise ratio and high 
spectral resolution data, as well as a better sampling, are re- 
quired to confirm this result. 

Third, some variability should also be expected due to out- 
moving clumps in the wind of O stars. In the case of CPup, 
spectral substructures similar to those observed in WR spectra 
are probably the consequence of excess emission from clumps 
(Eversberg et al. 1998). The variation across the line can reach 
an amplitude of about 5 % of the total line intensity. However, 
the substructures in ( Pup vary on time scales of a few hours 
only, which is too short to suit the variability time scales re- 
vealed by the study of our sample of Oef stars and is again at 
odds with the results of our intense monitoring of HD 14434. 
Finally, we should also consider the possibility of the propaga- 
tion of a clump (as discussed in the previous paragraph) across 
a putative confined structure within the wind (as discussed in 
the previous section) or a density wave. Such a scenario could 
lead to modulations in the profiles of the lines formed within 
the confined wind. The interaction of such perturbations with a 
corotating structure should be considered in future theoretical 
developments, in order to establish its impact on line profiles. 
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7. Conclusions and future work 

We have performed a spectral analysis of the Oef stars 
HD 192281, HD 14442 and HD 14434. The double peaked 
structure of the emission component of the He II A 4686 line 
of these three objects is obvious, and points to a possible 
corotating structure covering a significant fraction of the 
emitting region of this line in the stellar wind (Conti & Leep 
11974 ). We determined spectral types (respectively 05(ef), 
05.5ef and 06(ef) for the three stars) in agreement with 
former classifications. Our radial velocity measurements did 
not reveal any evidence of binarity on time scales of several 
days, or from one year to the next. The peculiar velocity of 
HD 192281 does not support a putative runaway status for this 
star. 

A detailed line profile variability analysis was performed 
on the strongest lines in the blue spectra of these stars. We have 
shown that the He II A 4686 line of HD 192281 and HD 14442 
displays the most significant variability, with time scales 
which are consistently found from one epoch to the other for 
both stars. We find however that the structure of the profile 
variability changes with time. The 'periodic' modulations of 
the profile, which are confined in a narrow wavelength domain, 
switch forth and back between the blue and the red wings, but 
never leave the central absorption component of the line^. The 
variability level itself changes also depending on the data set. 
Assuming that the periodicity of this variability is ruled by the 
rotation of the star, we propose rotational periods of 1.5 ±0.1 
and 1.6 ± 0.1 d for HD 192281 and HD 14442 respectively. 
HD 14434 displays the most complex behaviour on a yet 
undefined time scale. However, this star does not display 
any variations on time scales of a few hours, arguing against 
a pulsational origin for the observed line profile variations. 
The H/3 line displays lower level variations. In the case of 
HD 14442, the variations of Hf3 are correlated with those 
observed in He ii A 4686. 

In the framework of the most recent theoretical results con- 
cerning the magnetohydrodynamics of magnetic field-aligned 
rotating winds ("ud-Doula l2002> . we tentatively propose a sce- 
nario to explain some aspects of the observed line profile vari- 
abiUty of HD 192281 and HD 14442. In this scenario, the in- 
teraction of the magnetic field with the wind plasma leads to 
the confinement of a significant part of the stellar wind in the 
vicinity of the magnetic equator plane. While several aspects 
of the observed variability behaviour could be explained by a 
corotating structure resulting from this interaction, we note that 
this scenario cannot explain the wavelength confinement of the 
variability pattern nor the shift of this pattern from the blue to 
the red from one year to the next, if our sampling is adequate. 
The wavelength confinement of the variability around zero ve- 
locity we describe in this study suggests that the structures em- 
bedded in the stellar wind have a limited spatial extent, and 
the change in the part of the line which suffers this variability 

^ Spectral modulations on a time scale compatible with the rotation 
period, and located in the central region of the He II A 4686 line, were 
also observed in the case of Pup (Eversberg et al. ll998> . 



points to the fact that these structures are not static. Moreover, 
like in the case of BD H-60° 2522 (Paper I), the variations suf- 
fered by the line profiles of HD 14434 can not be explained by 
this scenario. 

The description of the variability of Oef stars constitutes a chal- 
lenge for theoretical modelizations. Further magnetohydrody- 
namical models including the interplay between NRPs, mag- 
netic fields and rotation need to be developed in order to clarify 
the complex behaviour of the lines in the spectra of Oef stars. 
On the other hand, from the observational point of view, a sim- 
ilar study of (Pup and ACep should be envisaged to check 
whether the variability time scales and the variation patterns 
reported in the literature are stable over longer time scales. We 
note also that the study of these stars should benefit strongly 
from the use of a wide-band spectrograph, as well as a denser 
time sampling. 
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